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Abstract 

Electron cooling is a well-established method to improve the phase space quality of ion beams in storage rings. In 
the common rest frame of the ion and the electron beam the ion is subjected to a drag force and it experiences a loss 
or a gain of energy which eventually reduces the energy spread of the ion beam. A calculation of this process is 
complicated as the electron velocity distribution is anisotropic and the cooling process takes place in a magnetic field 
which guides the electrons. In this paper the cooling force is calculated in a model of binary collisions (BC) between 
ions and magnetized electrons, in which the Coulomb interaction is treated up to second order as a perturbation to 
the helical motion of the electrons. The calculations are done with the help of an improved BC theory which is 
uniformly valid for any strength of the magnetic field and where the second-order two-body forces are treated in the 
interaction in Fourier space without specifying the interaction potential. The cooling force is explicitly calculated 
for a regularized and screened potential which is both of finite range and less singular than the Coulomb interaction 
at the origin. Closed expressions are derived for monochromatic electron beams, which are folded with the velocity 
distributions of the electrons and ions. The resulting cooling force is evaluated for anisotropic Maxwell velocity 
distributions of the electrons and ions. 
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1. Introduction 

In most experiments with particle beams a high phase space density is desired. In electron cooling of ion beams 
^ this is achieved by mixing the ion beam with a comoving electron beam which has a very small longitudinal 
momentum spread. In the rest frame of the beams the cooling process may be viewed as the stopping of ions in 
an electron plasma iQ-Ht]. More recently electron cooling has also been used in traps for precision experiments 
like CPT-tests with antihydrogen or planned QED-tests with highly charged ions in HITRAP [8]. In these 

applications the presence of strong external magnetic fields constitutes a theoretical challenge |9], as its influence 
on the cooling which the magnetized electrons exert on the ions (antiprotons) is not so obvious as earlier models 
might suggest. In the dielectric theory (DT) the drag on the ion is due to the polarization it creates in its wake. This 
can be either calculated in linear response (LR) fldjll!] or numerically by a particle-in-cell (PIC) simulation of the 
underlying nonlinear Vlasov-Poisson equation 1 12, 13]. While the LR requires cut-offs to exclude hard collisions of 
close particles the collectivity of the excitation can be taken into account in both approaches. In the complementary 
binary collision (BC) approximation the drag force is accumulated from the velocity transfers in individual collisions. 
This has been calculated by scattering statistical ensembles of magnetized electrons from the ions in the classical 
trajectory Monte-Carlo method (CTMC) II 1 3l-l 1 811 . and by treating the Coulomb interaction as a perturbation to the 
helical motion of the electrons, see Refs. ['19L'25'] for details. The observed cooling force F(v,) on an individual ion is 
obtained by integrating with respect to the impact parameter and the electrons velocity distribution. The ion velocity 
V; is measured with respect to the center of that distribution. As in electron cooler the electrons are accelerated from 
the cathode, their velocity distribution is flattened longitudinally, but the spread does not vanish. As the cooling force 
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on slow ions and therefore the coohng process depends critically on the details of the velocity distribution, a treatment 
employing a realistic velocity distribution is desirable. 

The purpose of this paper is the application of a second order perturbative BC model for calculating the magnetized 
cooling force on a uniformly moving individual heavy ion as well as on a heavy ion beam. In previous approaches lfl9[ 
I20I1 three regimes are identified, depending on the relative size of the cyclotron radius, the distance of the closest 
approach, and the pitch of the helix. The present paper is a continuation of our earlier studies in Refs. 112 1142511 where 
the second-order e nergy tr ansfer for an electron-ion colli sion is calculated with the help of an improved BC treatment 
which unlike Ref. il9ll20ll is valid for any strength of the magnetic field. 

In Section |2] we introduce a perturbative binary collision formulation in terms of the binary force acting between 
the ion and a magnetized electron, and derive general expressions for the second-order (with respect to the interaction 
potential) cooling forces. These expressions involve all cyclotron harmonics of the electrons' helical motion, and are 
valid for any interaction potential and any strength of the magnetic field. 

In Section [3] we turn to the explicit calculation of the second order cooling force without any restriction on the 
magnetic field in case of a regularized and screened interaction potential which is both of finite range and less singular 
than the Coulomb interaction at the origin and which includes as limiting cases the Debye (i.e., screened) and Coulomb 
potentials |2^ 27]. In addition we calculate the magnetized cooling force averaged with respect to the electron and 
ion beams velocity distribution functions. The cooling force for monoenergetic electrons is folded with an anisotropic 
velocity distributions which is typical for electron cooling of ion beams in storage rings, where the velocity spread is 
much smaller longitudinal than transverse to the magnetic guiding field. A similar anisotropic distribution is used for 
averaging with respect to the ion velocity distribution. Also asymptotic expressions for large and small ion velocities 
and strong and vanishing magnetic fields are given. 

Numerical results on the cooling force are presented in Section |4] using parameters of the ESR storage ring at 
GSI [28-30]. In particular, we compare our approach with the experimental data, the simplified treatment derived in 
Refs. L2(Ii25i1 and the model of Parkhomchuk Ii31il . 



2. Theoretical model 

2.1. Binary collision (BC) formulation 

We consider two point charges with masses m, M and charges -e, Ze, respectively, moving in a homogeneous 
magnetic field B - Eb. We assume that the particles interact with the potential -Z^U(r) with ^ = e^/A-nso, where 
eo is the permittivity of the vacuum and r = ri - r2 is the relative coordinate of the colliding particles. For two 
isolated charged particles this interaction is given by the Coulomb potential, i.e. t/c(r) = I/'"- In plasma applications 
the infinite range of this potential is modified by the screening. Then U may be modeled by Uoir) - e^''^'^/r with 
a screening length A, given e.g. by the Debye screening length Ad, see, for example 132 1. The quantum uncertainty 
principle prevents particles (for Z > 0) from falling into the center of these potentials. In a classical picture this can be 
simulated by regularizing U{r) at the origin, taking for example U^ir) = (1 - e^''^^)e^''l^ jr, where A is usually related 
to the (thermal) de Broglie wavelength [26,, 27i] . Here, however, the use of this regularized interaction essentially 
represents an alternative implementation of the standard (lower) cutoff procedure needed to handle the hard collisions 
in a classical perturbative approach. Hence we consider as a given constant or as a function of the classical collision 
diameter (see Section|4]i. 

In the presence of an external magnetic field, the Lagrangian and the corresponding equations of particles motion 
cannot, in general, be separated into parts describing the relative motion and the motion of the center of mass (cm) 



02511 ■ However, in the case of heavy ions, i.e. M » m, the equations of motion can be simplified by treating the cm 
velocity Vcm as a constant and equal to the ion velocity v,, i.e. Vcm = v, = const. Then the equation of relative motion 
turns into 

v(f) + uc [v(f) X b] = -u, [V,- X b] - -^f (r(f)) , (1) 

m 

where y{t) - r(t) - \e{t) - v,- is the relative electron-ion velocity, -Z/^f (r(f)) (f = -dU/dr) is the force exerted by 
the ion on the electron, cjc = eB/mis the electron cyclotron frequency. 
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It is now useful to introduce the velocity correction through relations 6\(t) - V(,(f) - \eQ(t) — v(f) - Vo(f)> where 
Vfo(f) and vo(f) are the unperturbed electron and relative velocities, respectively, with vo(f) = ro(f) - VeoCO - v,. 



ro (0 = Ro + v,.f + a [u sin (cdJ) - [b x u] cos (tOct)] , (2) 

(5v(f) + cjc [5v(0 X b] = — (r(f)) (3) 
m 

and 6\(t) ^ at f — > -oo. In Eq. (|2]) u = (cos (p, sin tp) is the unit vector perpendicular to the magnetic field, the angle 
ip is the initial phase of the electron's helical motion, vvy and Ve± (with Ve± ^ 0) are the unperturbed components of 
the electron velocity parallel and perpendicular to b, respectively, - v^yb - v, is the relative velocity of the guiding 
center of the electrons, and a - vv±/wc is the cyclotron radius. In Eq. (|2}, the variables u and Rq are independent 
and are defined by the initial conditions. In Eq. (O r(f) = rdt) - \it is the ion-electron relative coordinate. We also 
introduce the variable s - Ro± = Ro - nr(n,- ■ Rq) which is the component of Ro perpendicular to the relative velocity 
vector V,- with = Vj /v^. From Eq. ^ we can see that s is the distance of closest approach between the ion and the 
guiding center of the electron's helical motion. 

We seek an approximate solution of Eq. ^ in which the interaction force between the ion and electrons is consid- 
ered as a perturbation. Thus we are looking for a solution of Eq. Q for the variables r and v in a perturbative manner 
r = ro + ri + v = Vo + Vi + where ro(0, vo(0 are the unperturbed ion-electron relative coordinate and velocity, 
respectively, r„{t), v„(f) (n - 1,2,...) are the nth order perturbations of r(f) and v(0, which are proportional to Z". 

The equation for the first-order velocity correction is obtained from Eq. Q replacing on the rhs the exact relative 
coordinate r(f) by ro(f) with the solutions Vi(f) - ri(f) and 

ri(0 = — {-bP||(f) + Re [b(b • P^t)) - Px(f) + ' [b x P^(f)]]} . (4) 
m 

Here we have introduced the following abbreviations 

P|,(0 - r b ■ f (ro(T)) (t - t) dT, PAO = — r f (ro(T)) - l]dT (5) 

J -oo J -OO 

and have assumed that all corrections vanish at f — > -oo. As will be shown in the next Section, Eqs. (|2]i and (|4]l 
completely determine the second order cooling force on the ion. 

2.2. Second order cooling forces 

We now consider the interaction process of an individual ion with a homogeneous electron beam described by a 
velocity distribution function fiye) and a density ng. We assume that the ion experiences independent binary collisions 
(BCs) with the electrons. The total cooling force acting on the ion is then obtained by multiplying the binary force 
Z^f(r(f)) by the element of the electron relative flux ngVrd^sdt (where s is the impact parameter introduced above 
which is perpendicular to the relative velocity v, ) and integrating with respect to time and folding with velocity 
distribution of the electrons. The result reads 



F(v,) = Z^2n, r d\ef(ye)Vr f dh r f{r{t))dt 

%J %J \J —CO 



(6) 



and is an exact relation for uncorrected BCs of the ion with electrons. We evaluate this expression within a systematic 
perturbative treatment considering only the second order force F2 with respect to the binary interaction since the av- 



eraged first order force Fi vanishes due to symmetry reasons ||2]J-|25|]. Within the second order perturbation treatment 
the cooling force can be represented as: 

F2 =Z^V J'iv,/(v,)v,- Jt/kf/(k)k [k ■ ri(f)] e^'^^'W (7) 

Here we have introduced the two-particle interaction potential U (r) and the binary force f (r) has been written using 
Fourier transformation in space. Furthermore the factor e'^'^^'^ in the Fourier transformed binary force has been 
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expanded in a perturbative manner as ^ gik<^o{')^i + /jj . ri(f)], where ro(f) and ri(f) are the unperturbed and 

the first order corrected relative coordinate, Eqs. (0 and (|4]i, respectively, e.i. 



f (r(0) = -i (/kt/(k)ke 



'I 



i £/kJ/(k)k[l +ik-ri(f)]e 



ik-ro(0 



(8) 



From Eq. O it is seen that the second order cooling force, F2, is proportional to Z^. 

Substituting Eqs. (|4]i and (|5]l into Eq. (|7]) and writing the binary force in expression (|5) in terms of Fourier trans- 
formed potential results in 



m 



f d\J(\,)v, f dh r t/k£/k';7(k)f/(k')k f e*'°Wt/f f e'^' ''"^^^ dr 

%J %J %J %J —00 %J — oo 



(9) 



x\go(t-T) + — sin (w,.(f - t)) [1 - cos (tocit - t))] !> , 

OJr a>c 



where ^0 = (k ■ b)(k' ■ b), gi - (k k') - (k ■ b)(k' ■ b), ^2 = (k ■ [k' x b]). The time-integral in Eq. Q can be perfomed 
using the Fourier series expansion of the exponential function with Eq. (|2]l (see, e.g., Ref. Ii33i1 ). This yields 



— J d\J(\,)vrJ dh J dkdk'U(k)U{k')ke"-^^^"'-^o 

00 

Yj e'^"^"'^^e-'"' J„ ik^a) 7„, {k'^a) 5 (^„(k) + ^,„(k')) 



(10) 



x<^-- 



^0 



1 



1 



'g2 

lOJr 



^,„+i(k')-/0 ^,„-i(k')-/0 

1 1 



4,(k') - /O ^„,+i (kO - /O ^,,,-1 (k') - /O 



Here 7„ are the Bessel functions of the nth order, tan 6 - k^jk^, k\\ - (k ■ b) and fc^ are the components of k parallel and 
transverse to b, respectively, ^„(k) - ncoc -n k ■ v^, and is the initial phase of the electron as defined in the previous 
Section. Note that expression ( fTOt involves all cyclotron harmonics. 

Next, we integrate with respect to the initial phase ip and impact parameter s. For that purpose we recall that 
the volume element d\e can be represented in cylindrical coordinates as dy^ - dve\\Vi,j_dvi,±d(p, where v^w and Vgj_ are 
the electron velocity components parallel and transverse to b, respectively. The s-integration is enabled by using the 
relation e'^^° - e"'"'^""e"'^'^, where = k ■ n,-, = k - n,(k ■ n, ), i.e. the component of k parallel and transverse to 
iir. Performing now the <p and s-integrations results in 



F, = - 



(2ny Z-'^'^n 
2m 



00 Jo 



fiVe\\,Ve±)Ve±dv^ 



f 



dk\U(ktk 



(11) 



^ jl(k^a)\kl6'(U^))+^[6((„M)-Hi.-i(k))] 



where the prime indicates the derivative with respect to the argument. For deriving Eq. (fTTt we assumed an axially 
symmetric velocity distribution f{\e) = five\\, Ve±) and used 6(k\i)6(kj^) = S(k). 

The n-summation in Eq. ( fTTT l can be done using the summation formula for 2„ e'"^7^(z), see Ref. ll33ll . We then 
obtain 



F, = 



(iTTTzYn 



f 

Jo 



■ J dVe\\ J f (VeW, Vei) Ve±dVe± J dk\U (k)|^ k 



(12) 



k^ + k^^ 



sin (ciJct) 



Jq pfej^fl sin —^j sin (k ■ v, f) tdt. 



This is a general expression for the magnetized coohng force acting on an individual ion. It has been derived within 
second order perturbation theory but without any restriction on the strength of the magnetic field B. 

Some limiting cases can be easily extracted from Eq. (fTZt . For instance, at vanishing magnetic field sin(wcO/('^cO 
1 and the argument of the Bessel function should be replaced by kj_Ve^t. In the presence of an infinitely strong mag- 
netic field, however, the term in Eq. (fTZb proportional to k]_ and the argument of the Bessel function vanish since the 
cyclotron radius a — > 0. 



3. Cooling force for a regularized and screened Coulomb potential 

In electron cooling of ion beams the velocity distribution of the electrons is anisotropic. It is usually modeled by 
a two-temperature anisotropic Maxwell distribution. As the ion beam had a radial position some mm off the electron 
beam axis in the experiments lisUsoll. an additional transverse cyclotron velocity of the electrons has to be taken 
into account, see 1.30,1 for details. The velocity distribution relevant for the averaging in Eq. (fT2t is thus given by 

where ^ Vth± in the present case and the thermal velocities are related to electron temperatures by v^^j^^ - T^/m, 
^rtiii ~ T\\/m (here the temperatures are measured in energy units). In this case the transversal (F^ = F - bFy) and 
longitudinal (Fy = b ■ F) components of the cooling force (we drop here the index 2 for simplicity), Eqs. (fT2l i and 
(fT3l l, after velocity integrations |33| can be represented in the forms 



(v,-) 



f 

Jo 



-'^kla^^^-kU-{i-cost)L2 I /-2si"nj ^iCos(V/||f)-/i(M,±r) 
' " ^ t i\ sin(fe||fl,||f)7o(^±fl/j.O 



with Fj^(v,) = — f'j.(v,). Here we have assumed a spherically symmetric potential C/(k) = U{k) and have introduced 
the thermal cyclotron radii of the electrons Oj^ - v'th^/wc, a\\ - I'thn/'J^c, and 0,-^ - Vj^/ciJc, a,|| = vm/cjc, - v,|| - Ve- 
in general the cooling force is thus anisotropic with respect to the ion velocity v,-. 

For the Coulomb interaction U{k) - Ucik), the full 2D integration over the s-space results in a logarithmic diver- 
gence of the k-integration in Eqs. (fTTl l and (fT2l l. To cure this, cutoff parameters fcmin and fc,nax must be introduced, see 



1 22l 12511 for details. Instead of doing so, we here employ the regularized screened potential U(r) = UR(r) introduced 



in Section im with the Fourier transformed 



W)--^! ,, \_o - ^ ,_J . (15) 



where d'^ - A'^ + 

Up to now we have considered the magnetized cooling force acting on the individual ion which interacts with 
an electron beam with anisotropic velocity distribution. However, in some experiments the measured longitudinal 
cooling force represents an average over the drag forces on individual ions. For a comparison with our theoretical 
model the cooling force is thus interpreted as the average (F||(v,)) = IF of the component F||(v,) of the drag force (fl4l 
parallel to the beam axis (and the magnetic field) over the ion distribution /i(v,||, v,x) in the beam (see, e.g., 1i34h361 ). 
that is, 

dviW I fiivi\\,Vi±)Fn{vin,Vi^)vi^dvij_. (16) 

CO Jo 

Here we model the ion beam by the anisotropic Maxwell distribution 

(27r) ' o-^a-ii 
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where cr]_ - (l/2){vf_J = Tij_/M, cr^ - (v^|) - v|| - Tm/M with the effective transversal (r,x) and longitudinal (7^) 
temperatures of ions. Here M is the ion mass and v^n is the average cm velocity of the ion beam with respect to the 
electron beam. 

This yields the averaged cooling force (fT6t by substituting Eqs. (fT4l i and STJl into Eq. ( fTSI l and then integrating 
over Vi± and v,|| as 



9(u) = — 



"thii 



Jo Jo 



(k) k^dk^ 



(18) 



with 



, _ at 
D (f) = (5 + T — sm — 

I at 2 



(19) 



The dimensionless parameters introduced in Eq. ([18) are defined as m = (v,||-Vc)/ V2vthj|, a - Wc/}/vth||, 5y - l+cr^/Vjy|, 
(5 = cr^ /vth||, and r = T_i/T^\ is the anisotropy parameter of the electron beam. 

Finally substituting the interaction potential (fTsT i into Eq. ( fTST i and performing the A;||-integration we arrive, after 
lengthly but straightforward calculations, at 



-J(m) 



4 VizV^'n, 



^h|| 



. n dt_ 

"Jo f Jo 



fl'^0(iA(f,f))exp 



(20) 





L 2«2^2x 


2^^ sin {at) 


[i^^ (^) 


I ^^(f)J 


^ Qit,0 at 



Here P{0 = (<5^f' + 1 - QiUO = £>(f)^' + 1 - 'A(f,0 = (fV2)(l - i^)!^, 

2 



O (z) = e--^ + e-^'-^ - —1- i (e-^ - g-^'-^) , 



(21) 

- 1 z ^ ' 

where >c - Ajd - \+ AjX. Equation (l20l l is the main result of this paper. We next consider some Umiting cases of this 
expression. 

3.1. High-velocity and B — > limit 

In the high-velocity limit assuming "v/n > (Wc/l, Vth||;±, cry-j^) only small t contribute to the cooling force (120) 
due to the short time response of the electrons to the moving fast ions. In this limit we have sin(a'f)/a'f 1 and 
Q{t, 6\i:^ + 1 - where 6\ - D(Q) = 6^ + t. The remaining f-integration can be performed explicitly (see 



Appendix A for details) and Eq. (l20l i turns into 



- J(m) ^ T A(>f) 



"thii 



Jo 



exp 



■2„2 



Here as in Refs. ll22l425[l we have also introduced the generalized Coulomb logarithm (see [Appendix A i 

>^2 + 1 



■ In X - 1 . 



(22) 



(23) 



Note that Eq. (l22l i does not depend on the magnetic field, oJc, as natural consequence of the short time response of the 
magnetized electrons and is also valid for vanishing magnetic field. 
A further increase of the ion beam velocity yields 



-J(m) ^ 



2nZ^^*n, 



A(^)^ 



erf 



(24) 



2nZYneAi>c) 



"thii 



i.e. the force decreases as ^F(m) ~ u with the beam velocity. Here erf(z) is the error function. 



3.2. Strong magnetic field 

In the limit of very strong magnetic fields with a>cA » (v/y, Vth||;±, cr||;±) assuming that sin(Q'f)/af 
6^^^ + \ - from Eq. ( l20l i after lengthly but straightforward calculations we obtain 



and g(f, 



ljn7?ln^ 
- J{u) T A (K) 



u I 

Jo 



[(52^2+52(l_^2)]3/2 



exp 



2/2 



die +5^ {I 'C-) 



In particular, for the high-velocity limit with cjcA » m » 1, Eq. dZST l becomes 



- J (m) A i x) — 

mv:.„ M 



erf 



2 M / 2m2 



(25) 



(26) 



'thii 



■A(>^)-. 



There is an important difference if we compare Eqs. (l25T l and (|26] | with Eqs. (l22l l and (l24l l. respectively. The force 
( l25b decays as 3^(m) ~ m"^ much faster than in Eq. ( l24l i. The velocity of the beam in Eq. (l26b is large but is restricted 
to the value Wf/i, i.e. 1 <K m <«; Wc/l. Thus it cannot be arbitrary large. The velocity in Eqs. (l22i and (l24l i is arbitrary 
large but now restricted below, v/y s> UcA, i.e. the magnetic field there cannot be arbitrary large. 

Considering on the other hand also the case of small velocities m <K 1 at strong magnetic fields, Eq. (l25T l becomes 



8 V7rZV*ne 



"thir 



where 



2(1 



+ 2- 



3jc 



p{x) 



arccos x 
ln(x + Vx^ - 



1), 



X < 1 

X > 1 



(27) 



(28) 



(29) 



As expected the low-velocity cooling force Eq. dZTl l strongly depends on the details of the distribution functions of 
electrons and ions. 

More generally, a magnetic field is considered as strong, if the cyclotron radius a is smaller than the screening 
length A of the effective interaction between the electrons and the ions. Then the drag is dominated by the contributions 
from the adiabatic collisions with a < A. As the distance of closest approach is inversely proportional to the square 
of the relative velocity between the guiding center of the electron and the ion the influence of the magnetic field via 
adiabatic collisions tends to become more and more important for v, — > 0. 



3.3. Adjustment of the effective interaction 

For application of our present results, which are obtained with the statical screened interaction U^(r), to the elec- 
tron cooling problem the static screening length A has to be replaced by an appropriately chosen velocity-dependent, 
dynamical one. For a known average temperature T - \{T\\ + IT^} of the electrons with average thermal velocity 
^'th = (T/my^- this canbedone, for instance, by choosing the dynamical screening length /i(v/||) = /iD[l +(v,||/vth)^]'^^, 
where is the averaged Debye radius - Vth/(^p and aip is the electron plasma frequency. 

Next we have to specify the cutoff parameter A which is a measure of the softening of the interaction potential 
at short distances. As we discussed in the previous section the regularization of the potential (fTsT i guarantees the 
existence of the s-integrations, but there remains the problem of treating hard collisions as accurate as possible. For a 
perturbation treatment the change in relative velocity of the particles must be small compared to and this condition 
is increasingly difficult to fulfill in the regime — > 0. This suggests to soften the potential near the origin the 
more the smaller Vr is. In fact the parameter A should be related to the de Broglie wavelength which is inversely 
proportional to v^. Here, within a classical picture, we employ for the perturbative treatment a dynamical cutoff 



parameter >t(viii} = I + A(vi\i) / A(vi\i) ||23l|24[], where A^Cviw) = Cbl(vii\) + Al with bo(yi\y) = \Z\f /t7i(vf.. + \ 



+ Vjjj), where 



10' 




relative ion velocity (m/s) relative ion velocity (m/s) 



Figure 1 : Longitudinal cooling force (in eV/m) for C^* (left panel) and U'^* (right panel) fully stripped ions as function of the relative ion velocity 
(in m/s) with respect to the rest frame of the electron beam. The theoretical cooling force l l2Qt is calculated for Aq = 10"' m and for an electron 
beam with = 10^ cm"^, Tx = 0. 11 eV and = 0.1 meV in a magnetic field of B = (solid line), 0.02 T (dashed line), 1 T (dotted line), B = oo 
(dash-dotted line). The ion beam is characterized by the distribution fry = 0, o-j_ = 3<v,x), Si = 0.2 mrad (see the text for details). 



/to is some constant cutoff parameter, and bo is the averaged distance of closest approach of two charged particles in 
the absence of a magnetic field. Also in /t(y,||) we have introduced a fitting parameter C - 0.292. In Refs. l23l 12411 
this parameter is deduced from the comparison of the second order scattering cross sections with an exact asymptotic 
expression derived in Ref. iIjtIi for the Yukawa-type (i.e., with /t — » 0) interaction potential. As has been shown 
in Refs. ||23, 24] the second-order cross sections for electron-electron and electron-ion (with and without cyclotron 
gyration of the ion) collisions with dynamical cutoff parameter /t(v,||) excellently agrees with CTMC simulations at 
high velocities. As a consequence, the generalized Coulomb logarithm A(>c) determined by Eq. (l23l l depends now 
on the ion beam velocity viy and behaves at high-velocities as A(k) ^ Inn - I ^ ln(v,||/(yplo) - 1- This velocity 
dependent behavior of A(k) must be taken into account when considering the asymptotic expressions ( l22l i and (l25T l. 



4. Results and discussion 

With the theoretical formalism presented so far, we now study the cooling forces on the ions resulting from our 
analytical approach, Eq. (|20] |. and compare them with available experimental data. But before starting this, we first 
consider some general properties of Eq. (l20i . in particular, the effect of the magnetic field and of a variation of the 
ion and electron distributions on the cooling force. In Fig. [T]the forces are plotted vs ion beam relative velocity for 
different values of the magnetic field. The two limiting cases of vanishing (B - 0, solid lines) and infinitely strong 
(B = oo, dash-dotted lines) magnetic fields are obtained from Eqs. (l22l i and ( l25T l, respectively. The density and the 



temperatures of the electron beam are the same as in the experiments at ESR storage ring 11281 - 13011 (see below for 
details). The parameters cry and o-± of the ion beam are given in units of the quantity (v,x) introduced below and are 
typical for many electron cooling experiments (see, e.g., Refs. UBS). It is seen that the magnetic field increases 



the cooling force 3^(u) at low velocities while reducing it at high-velocities. As expected the cooling force is more 
sensitive to the magnetic field at small u where a rather weak field with B - 0.02 T produces some deviations from 
the B = regime. At high-velocities and at strong magnetic field (B = 1 T) the cooling force strongly deviates from 
the extreme case with B - oo, which is, however, not accessible for the present experiments at storage rings. Also 
we would like to note that at very high-velocities the cooling force ( l20l i systematically deviates from Eq. (l22t which 
is the leading order term 0(yT^) of high-velocity expansion of Eq. ( l20l i. This is because at m » 1 the next term 
with 0(v^|'*) in high-velocity expansion of Eq. ( |20] | is proportional to (wc/lvth||/v^|)^ and at high-velocities behaves as 
(Wc/wp)^(v'th||/vi||)^ ~ 0(y'^^) due to the dynamical screening introduced above. Thus, at m » 1 and at finite magnetic 
field Eq. (l20t deviates systematically from the regime of vanishing magnetic field, Eq. ( l22l i. 




relative ion velocity (m/s) relative ion velocity (m/s) 



Figure 2: Same as in Fig.^for C^* ions and for fixed magnetic field field with B = 0.1 T varying the shape of the distribution function of ions (left 
panel) and electrons (right panel). Left panel, (Ty = 0, the lines without symbols and o"|| = 2{v,x), the lines with symbols, ctx = 0.1{v,x) (solid 
line), <T± = (v,x> (dashed line), ctx = 3.5(v,x) (dotted line). Right panel, iT|| = 0, <t± = 3.5(v,x), = 0.1 meV, the lines without symbols and 
= 1 eV, the lines with symbols. T± = OA meV (solid line), T± =0.1 eV (dashed line), T± = I eW (dotted line). 



Next, in Fig.|2] the cooling forces on C^^ ions are plotted vs ion beam relative velocity for different shapes of the 
distribution functions of ions (left panel) and electrons (right panel). As expected at high-velocities the forces do not 
depend on the details of the distribution functions of the electron or ion beam. In the low-velocity regime 3^{u) is 
sensitive to both distributions showing, however, essential sensitivity to the transversal velocity spread of ions at small 
cr||, i.e., when the ion beam is nearly monochromatic in longitudinal direction (Fig.|2] left panel). As a general property 
of expression (l20l i we also mention the broadening of the profile of the cooling force at small cr|| with increasing the 
velocity spread o-± in transversal direction. Similarly, for a given ion beam (i.e. for a fixed cry-x), the force profile is 
broadened with decreasing longitudinal energy spread Ty of the electrons. 

Measurements of the cooling forces have been performed at several storage rings, like e.g. at the ESR at GSI 1128 - 
I30I1 . In these experiments a so-called cooling force is extracted, which can be viewed as a stopping force averaged 
over the ion distribution in the beam and the electron distribution. As an example we focus on the measurements 
of longitudinal cooling forces for diff'erent fully stripped heavy ions as conducted at the electron cooler of the ESR 
storage ring. Two different methods have been used here to determine the cooling force. At low ion velocities the 
cooling force is extracted from the equilibrium between cooling and longitudinal heating with rf noise. At high relative 
velocities between the rest frames of the beams the cooling force is deduced from the momentum drift of the ion beam 
after a rapid change of the electron energy. Details of these methods are given in ll28i I29I1 . The measured cooling 
forces are shown in Figs.|3]and|4](filled circles). 

The electron beam in these experiments has a density of - 10^ cm"^ and can be described by an anisotropic 
velocity distribution (fT3l l with T± - mv^^^ ^ 0.11 eV and T\\ - wv^y - 0.1 meV. As the ion beam had a radial 
position some mm off" the electron beam axis injhese experiments, an additional transverse cyclotron velocity of 
the electrons has to be taken into account, see [30] for details. The velocity distribution relevant for the averaging in 
Eq. ( fT2b is thus given by Eq. ( fT3T l. where - Vth± in the present case. The strength of the magnetic guiding field 
was B = 0.1 T. The measured longitudinal cooling force represents an average over the stopping forces on individual 
ions. For a comparison with the theoretical model (l20t the cooling force is thus interpreted as the average {F\\) of the 
component F\\ of the stopping force (fT2l i parallel to the beam axis (and the magnetic field) over the ion distribution 
fi{vi\\, v,x) in the beam (see also [34-36]). 

For low ion velocities this average is taken with respect to the transverse ion velocity only and the cooling force 
depends on the parallel ion velocity, i.e. {F\\) - (Fy)(v,||). In the experimental procedure used for high ion velocities 
the cooling force is an average over the complete ion distribution. This average {F\\) - (f'y)((v,y)) depends now on 
the velocity of the cm of the ion beam relative to the rest frame of the electron beam (v,||). Both velocities are denoted 
as relative ion velocity in Figs. |3]and|4] To perform the average the distribution fiiv^, v,x) must be known. However, 
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Figure 3: Longitudinal cooling force (in eV/m) for various fully stripped ions as function of the relative ion velocity (in m/s). Filled circles: 
experimental data from measurements at the electron cooler of the ESR storage ring I28l430ll. Dotted curves: binary collision approximation 
proposed in (l^l . Dashed curves: empiric formula 1 1301 for the cooling force as proposed in l3lll . Solid curves: equation ilOi . The theoretical 
descriptions of the cooling force are calculated for an electron beam with = lO'' cm"^, T± = 0.11 eV and = 0.1 meV in a magnetic field of 
B = 0. 1 T, and are fitted to the experimental results at low relative velocities by treating the transverse ion velocity (v;^) (for dotted and dashed 
curves) or the quantities (Ty-x (for solid curves) as a free parameter (see the text for details). 



in Refs. 1I28H30I1 this distribution was not determined in detail, but there exists an estimate of the beam divergence 
(0,) < 0.5 mrad 1E9I1 . This yields after transformation to the rest frame of the ion beam for the transverse ion velocities 
{vi±} - fiyc{0i), where /?, y are the relativistic factors related to the beam velocity in the lab frame and c is the speed 
of light. For the measurements at hand with an ion energy of 250 MeV/u (fi - 0.615, y - 1.268) this results in 
<v',±) < 10^ m/s. 

In our previous simplified model ['25'] due to the lack of a more detailed knowledge about the ion distribution, the 
average with respect to /;(v,||, v;±) is replaced by {F\\) - F\\{\'n_ - (v,j.), v',||), where v,j. = (v,j_) was treated as a free 
parameter to fit the BC stopping force to the experimental data. As F\\(vn_, v,||) is rather sensitive to a variation of v,x 
at low parallel velocities v,|| this fit is done for the linear increase of the cooling force at low relative velocities. The 
related values for v/j^ = (v,j^) are in the range v,x/vth|| = 10-17 which corresponds to (0,) ^ 0.2 - 0.3 mrad in good 
agreement with the estimated beam divergence. The resulting theoretical predictions are given by the dotted curves 
in Figs. [3] and m They agree well with the experimental data at low and high velocities but overestimate the cooling 
force at medium velocities Vthy - 4.2 x lO"' m/s < v,|| < Vth± = 1.4 x 10^ m/s. 
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Figure 4: Same as in Fig.|3]but for U'^^ ions. 



The experimental data are also compared to an empirical formula for the cooling force 

\nnj?ii'' v,|| / <,w^) \ 
F\\(yi) = —I 5 5—777 In 1 + — (30) 



as proposed by Parkhomchuk Bill . Here (imin) = jmiy^ + and (imax) = (i'^ + ^rti)'^^/'^p minimal and 

maximal impact parameters, - v±^/a)c is the cyclotron radius of the electrons, and Vgff is an effective electron 
velocity related to the transverse magnetic and electric fields in the electron cooler, see pll], which can be viewed 
as a fitting parameter The force (l30t , which is the stopping force on a single ion, must also be averaged over the 
ion distribution like in ( fTSI l. As above the force (l30t with an average (v,j^) has been used instead. If the additional 
parameter is chosen as Veft = Vth|| (that is, rather small) the experimental data at low velocities are fitted by nearly 
the same values for the beam divergence (0,) as before. The resulting cooling force is shown by the dashed curves 
in Figs. |3] and |4] The agreement with the experimentally deduced cooling force is as good as for the simplified BC 
force proposed in Ref. [25] (dotted curves) at low and high velocities. At medium velocities the deviation from the 
experimental results is even larger, in particular for v,|| ^ Vth±, compared to the simplified but more detailed BC 
treatment of Ref. |25]. 

Now we turn to the present expression for the cooling force (EOl i which is shown as solid curves in Figs. |3] and |4] 
assuming the same beam divergence as in previous cases (i.e. (0,) ^ 0.2 mrad). The velocity spread of the ion beam 
in transversal direction used in obtaining the solid curves is 3.5/5yc{6i) ;$ < A.5[iyc{9i) while the spread in the 
longitudinal direction is typically cr|| < IO^^ctj^ as it usually occurs in many experimental situations (see, e.g.. 



and references therein), in particular at ESR storage ring [28-30]. It is seen that the agreement of the more detailed 
Eq. ( I2OI 1 with the experimental cooling force is essentially improved in the whole relative velocity range compared 
to the simple models considered above (dotted and dashed curves). We mainly ascribe the deviations of the present 
model (solid curves) from the ESR data to the rather unknown distribution function of the ions in the beam which 
has been modeled here in the form of an anisotropic Maxwell distribution (fTTT i. Indeed the actual velocity spread in 
ion beams may essentially differ from the Maxwellian (flTt and, in particular, in some cases the recorded profiles are 
parabolic rather than Maxwellian 1,34-36.1 (see also Ref. |3]). For a comprehensive comparison with the measurements 
and a critical evaluation of the theoretical approaches a detailed knowledge of the ion distribution is indispensable. 



5. Summary 

In this paper we have calculated the cooling force on ions in a model of binary collisions (BC) between ions and 
magnetized electrons within second order perturbative treatment. The calculations have been done with the help of an 
improved BC theory which is valid for any strength of the magnetic field. The cooling force is explicitly calculated 
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for a regularized and screened Coulomb potential. Closed expressions have been derived for monochromatic electron 
beams, which have been folded with the velocity distributions of the electrons and ions. The resulting cooling force 
is evaluated for anisotropic Maxwell velocity distributions of the electrons and ions. In addition, a number of limiting 
and asymptotic regimes of low- and high-velocities as well as vanishing and strong magnetic fields have been studied. 
The given results show that the present model of the cooling force is very sensitive to the velocity spreads of the 
electrons and ions at small relative velocities. 

We have also compared our BC model with the previous theoretical approaches of Refs. 11251 bill as well as with the 



experiments performed at the ESR at GSI ll28[ - l30ll . As demonstrated, the agreement of Eq. ( l20l l with the experimental 
cooling forces considered over the whole relative velocity range is much better in comparison to the simplified BC 
model derived in Ref. li25ll and Eq. (|30] | proposed by Parkhomchuk. The deviations of Eq. (l20b from the ESR data, 
which can be seen in Figs. [3] and H) are ascribed to the deviations of the model distribution function (fTTT i from the 
experimental distribution of the ion beam which is not precisely known. 

As the main goal of this paper we suggest a more advanced analytical model for calculations of the cooling force 
which is appropriate for modeling many experimental situations with moderate or strong magnetic guiding fields. 
In principle, the average involved in Eq. (fTST i can also be done numerically with recorded ion beam distributions or 



analytically using other ion distributions like e.g. the parabolic distribution function as it occurs in CELSIUS Il35[ |3 



Systematic comparisons for different distribution functions and other experiments on electron cooling as well as with 
Monte Carlo numerical simulations are currently in progress and the results will be reported elsewhere. 

Finally we would like to mention that our results for the cooling forces F||(v;) and 3^(u) can be tabulated in a 
suitable manner to be used as input for simulations of electron cooling using the BETACOOL package issjlsoll . 
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Appendix A. Derivation of the Coulomb logarithm 

In this Appendix we evaluate briefly the f-integral remaining in Eq. ([20} in the limit of high-velocity of the ion 
beam. This integral is given by 

f" dt 1 r°° dz 

-0(<A(f,^))= - (A.l) 
Jo t 2 z 

where e ^ +0 and the function (i>(z) is determined by Eq. (ISTT i. The second relation in Eq. (lA.H is obtained by 
introducing the new variable of integration (f^/2)(l/^^ - 1) = z. Note that <i>{z) ~ at z — > and the second integral 
in Eq. (lA.ll i is convergent at vanishing z- However, we have introduced an infinitesimal e — > +0 as the lower limit 
of integration which allows to perform separately z-integration of each term of the function <l>(z). The calculation is 
straightforward and the final result reads 

"^'^^ -[£i(e)-£i(>^2e)]-l, (A.2) 



2(>*2_ 1)^ 



where £'i(z) is the exponential integral. At small argument (z — > 0) this integral behaves as £'i(z) - ln(l/z) - y 113 311 . 
where y is Euler's constant. For e — > +0, Eq. ( IA.2l i thus turns into the generalized Coulomb logarithm A(>c) given by 
Eq. dH. 
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